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ABSTRACT 

We report the results of an analysis of all Spitzer/MIPS 24 /xm observations of HD 209458b, one 
of the touchstone objects in the study of irradiated giant planet atmospheres. Altogether we analyze 
two and a half transits, three eclipses, and a 58-hour near-continuous observation designed to detect 
the planet's thermal phase curve. The results of our analysis are: (1) A mean transit depth of 
1.484% ± 0.033%, consistent with previous measurements and showing no evidence of variability in 
transit depth at the 3% level. (2) A mean eclipse depth of 0.338% ± 0.026%, somewhat higher than 
that previously reported for this system; this new value brings observations into better agreement 
with models. From this eclipse depth we estimate an average dayside brightness temperature of 
1320 K ± 80 K; the dayside flux shows no evidence of variability at the 12% level. (3) Eclipses in the 
system occur 32 ± 129 s earlier than would be expected from a circular orbit, which constrains the 
orbital quantity ecosw to be 0.00004 ± 0.00033. This result is fully consistent with a circular orbit 
and sets an upper limit of 140 m s^^ (3(7) on any eccentricity-induced velocity offset during transit. 
The phase curve observations (including one of the transits) exhibits an anomalous trend similar to 
the detector ramp seen in previous Spitzer/IRAC observations; by modeling this ramp we recover 
the system parameters for this transit. The long-duration photometry which follows the ramp and 
transit exhibits a gradual ~ 0.2% decrease in flux over ^ 30 hr. This effect is similar to that seen in 
pre-launch calibration data taken with the 24 /im array and is better fit by an instrumental model than 
a model invoking planetary emission. The large uncertainties associated with this poorly-understood, 
likely instrumental effect prevent us from usefully constraining the planet's thermal phase curve. Our 
observations highlight the need for a thorough understanding of detector-related instrumental effects 
on long time scales when making the high-precision mid-infrared measurements planned for future 
missions such as EChO, SPICA, and JWST. 

Subject headings: transits — eclipses — infrared: planetary systems — planets and satellites: individ- 
ual (HD 209458b) — planetary systems — techniques: photometric — stars: indi- 
vidual (HD 209458b) 



1. INTRODUCTION 

Most known extrasolar planets were discovered via the 
radial velocity technique - in which the Doppler wobble 
of a star indicates an orbiting planet - and/or by the 
transit method - in which periodic dimming of a star in- 
dicates a planet that crosses in front of the stellar disk. 
Owing to the observational biases of these techniques, the 
first planets thus discovered were the large, massive ob- 
jects on few-day orbits commonly known as hot Jupiters 
(Mayor & Qucloz 1995; Henry ct al. 2000; Charbonneau 
et al. 2000). Their large sizes and high temperatures 
make these objects excellent candidates for the study of 
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their dayside emission when the planet is occulted by the 
star (Doming et al. 2005; Charbonneau et al. 2005), of 
their longitudinally-averaged global emission (Harring- 
ton et al. 2006; Cowan et al. 2007; Knutson et al. 2007), 
and of their atmospheric opacity via the wavelength- 
dependent flux diminution during transit (Seager & Sas- 
sclov 2000; Charbonneau ct al. 2002). These observations 
have led to measurements of atmospheric abundances of 
key molecular species (Madhusudhan ct al. 2011a), pos- 
sible non-equilibrium chemistry (Stevenson ct al. 2010), 
high-altitude hazes (Sing et al. 2009), and atmospheric 
circulation (Cowan & Agol 2011b). 

Any discussion of hot Jupiter atmospheres must nec- 
essarily mention two systems in particular. One, 
HD 189733, is the brightest star known to host a hot 
Jupiter (Bouchy et al. 2005). The other is HD 209458, 
the first known transiting planet (Charbonneau et al. 
2000; Henry et al. 2000) and the focus of this study. 
These are the two touchstone objects in the study of irra- 
diated giant exoplanets, both because they were discov- 
ered relatively early on and because they orbit especially 
bright (as seen from Earth) host stars. This last point in 
particular allows for especially precise characterization 
of these planets' atmospheres and permits observations 
which would provide unacceptably low signal to noise 
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ratios for fainter systems. 

1.1. The HD 209458 system 

The star HD 209458 is an F8 star roughly 15 % more 
massive than the Sun (Mazch ct ah 2000; Brown ct al. 
2001; Baincs et ah 2008), with an equivalent metallic- 
ity and slightly higher temperature (Schuler et al. 2011). 
It is orbited by HD 209458b, a roughly 1.4 0.7 Mj 
planet in a 3.5-day, near-circular orbit (Southworth 2008; 
Torres et al. 2008). The planet's parameters have been 
substantially improved upon since its initial discovery 
(Charbonneau ct al. 2000; Henry et al. 2000; Mazeh et al. 
2000). Two sets of more recent values (Torres et al. 2008; 
Southworth 2008) do not differ significantly, and we use 
the former's system parameters in our analysis when not 
making our own measurements. 

Infrared photometry during eclipses of HD 209458b 
measured from the ground (Richardson et al. 2003) and 
with Spitzer (Doming et al. 2005; Knutson et al. 2008) 
determines the planet's intrinsic emission spectrum, and 
is best fit by atmospheric models in which the planet's 
atmospheric temperature increases above ~ 0.1 — 1 bar 
(Burrows et al. 2007, 2008; Fortney et al. 2008; Mad- 
liusudhan & Seager 2010). Such temperature inversions 
are common on hot Jupiters, and a popular explanation 
requires the presence of a high-altitude absorber (e.g., 
Fortney et al. 2008; Burrows et al. 2008). The nature of 
any such absorber is currently unknown and the subject 
remains a topic of active research (Desert et al. 2008; 
Spiegel et al. 2009; Knutson ct al. 2010; Madhusudhan 
et al 2011b). 

If present, a high-altitude optical absorber is expected 
to absorb the incident stellar flux high in the atmosphere 
where radiative timescales are short and advection is 
inefficient (Cowan & Agol 2011a). Consequently, such 
planets are expected to exhibit large day/night temper- 
ature contrasts and low global energy redistribution de- 
spite circulation models' ubiquitous predictions of large- 
scale superrotating jets on these planets (Showman & 
Guillot 2002; Cooper & Showman 2005; Cho et al. 2008; 
Rauschcr et al. 2008; Showman ct al. 2009; Dobbs-Dixon 
ct al. 2010; Burrows ct al. 2010; Rauschcr & Menou 2010; 
Thrastarson & Y-K. Cho 2010; Hcng ct al. 2011b,a). 
Spitzer/IRAC observations of HD 209458b at 8 fj,m place 
an upper limit on the planet's thermal phase variation 
of 0.0022 (3cr; Cowan et al. 2007). Given the planet's 
demonstrably low albedo (Rowe et al. 2008) this limit 
is substantially lower expected if the planet has a low 
recirculation efficiency. In hot Jupiter atmospheres the 
dominant 24 /im molecular opacity source is expected to 
be H2O, but there is some tension between models and 
past observations at this wavelength (cf. Madhusudhan 
& Seager 2010). Thus our understanding of these plan- 
ets' atmospheres remains incomplete. 

Recent spectroscopic observations of HD 209458b dur- 
ing transit show a hint of a systematic velocity offset 
(2 ± 1km s^^) of planetary CO lines during planetary 
transit (Snellen et al. 2010). If confirmed, this offset 
would be diagnostic of high-altitude winds averaged over 
the planet's day/night terminator, and similar measure- 
ments at higher precision could one day hope to spa- 
tially resolve terminator circulation patterns and con- 
strain atmospheric drag properties (Rauschcr & Menou 
2012). However, small orbital eccentricities (specifically. 



nonzero ecosw, where uj is the longitude of periastron) 
can also induce a velocity offset in a planetary trans- 
mission spectrum (Montalto et al. 2011). It is thus 
convenient that precise timing of planetary transits and 
eclipses directly constrains ecoscj (Seager 2011, chapter 
by J. Winn). This provides a further motivation for our 
work: to more tightly constrain HD 209458b's orbit via 
a homogeneous analysis of a single, comprehensive data 
set. 

In this paper we analyze the full complement of data 
for the HD 209458 system taken with the MIPS 24 ^m 
camera (which we hereafter refer to simply as MIPS; 
Ricke ct al. 2004) on the Spitzer Space Telescope. MIPS 
has taken previous 24 /zm observations of exoplanetary 
transits (Richardson et al. 2006; Knutson et al. 2009a), 
eclipses (Doming ct al. 2005; Charbonneau et al. 2008; 
Knutson et al. 2008, 2009b; Stevenson et al. 2010), and 
thermal phase curves (Harrington et al. 2006; Knutson 
et al. 2009b; Crossfield et al. 2010). MIPS operations de- 
pended on cryogenic temperatures; since Spitzer's com- 
plement of cryogen has been exhausted there may be no 
further exoplanet measurements at wavelengths > 10 /im 
until the eventual launch of missions such as EChO, 
SPICA, or the James Webb Space Telescope (JWST). 
Our work here describes some of the last unpublished 
24 /im exoplanet observations, and a further motivation 
for our work is to inform the calibration, reduction, 
and observational methodologies of future missions' mid- 
infrared (MIR) observations. 

1.2. Outline 

This report is organized as follows: in Section 2 we de- 
scribe the MIPS observations and our approach to mea- 
suring precise system photometry. In Section 3 we de- 
scribe our efforts to understand the origin of instrumen- 
tal sensitivity variations apparent in the long-duration 
phase curve observations; these effects ultimately prevent 
any measurement of HD 209458b's thermal phase curve. 
However, we are able to recover the parameters of the ob- 
served transits and eclipses, and we present these results 
in Sec. 4 and 5, respectively. Combining the results of 
these two analyses allows us to constrain the planet's or- 
bit (i.e., ecosw), and we discuss the implications of this, 
and of the total system flux, in Section 6. We summa- 
rize our conclusions and present some thoughts for future 
high-precision MIR observations in Section 7. 

2. OBSERVATIONS AND ANALYSIS 

2.1. Observations 

We reanalyzed all observations of the HD 209458 sys- 
tem taken with Spitzer's MIPS 24 /im channel: analysis 
of one transit, two eclipses, and the long-duration phase 
curve observations have remained unpublished until now. 
Altogether, we used the data from Spitzer Program IDs 
3405 (PI Seager; published in Deming ct al. 2005), 20605 
(PI Harrington; published in Richardson et al. 2006), and 
40280 (PI Knutson). Table 1 lists the observatory pa- 
rameters used for each set of observations. Collectively 
these data comprise 2.5 transits, three eclipses, and a 
58-hour set of near-continuous observations designed to 
detect the planet's thermal phase curve. 

2.2. Data Reduction 
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Unless stated otherwise we use the same methodology 
to reduce our data as described in (Crossfield et al. 2010, 
hereafter CIO), performing PSF-fitting photometry us- 
ing a 100 X super-sampled MIPS PSF^ modeled using a 
6070 K blackbody spectrum simulated at the center of 
the MIPS field of view. We vary the size of the synthetic 
aperture used to calculate our PSF-fitting photometry, 
and find that a square, 21x21 pixel aperture minimizes 
photometric variations. During MIPS observations the 
target star is dithered between fourteen positions on the 
detector (Colbert, J. 2011, Section 8.2.1.2), and we fit 
the data from all dither positions simultaneously as de- 
scribed below. 

As noted by CIO, the MIPS 24 /^m detector appears 
to suffer from low-amplitude temporal variations in the 
diffuse background, presumably owing to small amounts 
of scattered light in the instrument. Because this could 
affect the flat-fielding performed by the MIPS reduction 
pipeline, we create an empirical flat field by taking a 
pixel-by-pixel median of all the individual frames after 
masking the region containing the target star. After con- 
structing this flat field we extract photometry (a) after 
subtracting the master flat field from each frame, and 
(b) after dividing each frame by the normalized-to-unity 
master flat field. Both of these give photometry that is 
very slightly less noisy (RMS reduced by < 1%) than 
photometry that does not use an additional flat fleld 
correction. Subtracting by the empirical flat-field prior 
to computing PSF-fitting photometry results in a lower 
residual RMS and so we use this approach for all our 
data; ultimately our choice of flat field does not change 
our final results. 

We extract the heliocentric Julian Date (HJD) from 
the timing tags in each BCD data file, and then con- 
vert the HJD values to BJDtds using the IDL rou- 
tine hjd2bjd^° (Eastman et al. 2010). These new time 
stamps have an estimated accuracy of one second (East- 
man et al. 2010), which is small compared to our final 
ephemeris uncertainties of roughly one and four minutes 
for transits and eclipses, respectively. 

2.3. Approach to Model Fitting 

The MIPS dither pattern introduces systematic offsets 
of < 1% (Doming et al. 2005) in the photometry at each 
dither position. We follow the methodology of CIO and 
explicitly fit for this effect by multiplying the modelled 
flux for each visit at dither position i by the factor (l-fc^). 
We further impose the constraint that these corrections 
do not change the absolute flux level, and so define cq 
such that the quantity Hi (1 + '^0 equal to unity. We 
ultimately find that the Ci are similar, but not constant, 
from one epoch to the next. 

In all cases we determine best-fit model param- 
eters using the Python simplex minimization rou- 
tine scipy. optimize. fmin^^. We assess parame- 
ter uncertainties using a Markov Chain Monte Carlo 
implementation of the Metropolis-Hastings algorithm 
(analysis . genericjncmc^^), then take as uncertainties 

® Generated using Tiny Tim; available at http ; //ssc . spitzer . 
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Fig. 1. — MIPS 24 fim photometry of the HD 209458 system, 
showing the detector ramp (0-10 h), transit (5 h), and eclipse 
(48 h). For plotting purposes the data have been binned to lower 
temporal resolution. A slight 0.2%) flux decrease is apparent 
from 10-58 h. This could be influenced by planetary phase vari- 
ations, but the similarity to the purely instrumental effects seen 
in Figure 2 precludes an unambiguous distinction between the two 
effects. 

the range of values (centered on the best-fit value) that 
enclose 68.3% of the posterior distribution. We verify 
by eye that the Markov chains are well-mixed; the re- 
sulting one-dimensional posterior distributions are uni- 
modal, symmetric, and approximately Gaussian unless 
stated otherwise. 

3. CALIBRATION AND INSTRUMENT STABILITY 
3.1. The Ramp 

Before we present the results of our model fits, we dis- 
cuss two photometric variations that we conclude to be 
of instrumental origin. The HD 209458 system fiux mea- 
sured from our 2008 observations, shown in Figure 1, 
exhibits a steep rise during the first 10-12 hours in which 
the measured system fiux increases by ~ 2%. This ramp 
appears similar to that seen in photometric observations 
taken with Spitzer/IRAC and Spitzer/IRS (Cliarbon- 
ncau ct al. 2005; Doming ct al. 2006; Knutson et al. 
2007). The IRAC ramp is the better studied, and is 
thought to result from charge-trapping in the detector 
(cf. Knutson et al. 2007; Agol et al. 2010). According 
to this explanation, a substantial fraction of photoelec- 
trons liberated early in the observations become trapped 
by detector impurities, resulting in a lower effective gain 
for the detector. Eventually all charge-trapping sites be- 
come populated and the detector response asymptotes 
to a constant level. As the IRAC 8 /im, IRS 16 ^m, and 
MIPS detectors are all constructed of Si:As it is conceiv- 
able that the MIPS ramp we observe has a similar origin 
in charge-trapping. 

To test this hypothesis, we look for evidence of per- 
sistence in our data. Using all frames taken at the 
second dither position we compute the median image 
from each of several Astronomical Observing Requests 
(AORs). An AOR is a Spitzer logistical unit compris- 
ing some dozens of frames; in our data set each AOR 
lasts approximately 3 hr. We see faint afterimages at the 
other thirteen dither positions when we subtract the first 
median AOR frame from the final median AOR frame 
(taken ^ 56 hr later; cf. Figure 1), which suggests that 

python/ 
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the level of persistence (a byproduct of charge trapping) 
increases over the course of the observations. These af- 
terimages are much fainter when comparing data from 
the first and second AORs (separated by 3.6 hr), con- 
sistent with the conclusion that the level of persistence 
does not saturate to a constant value on these short time 
scales. The afterimages are not apparent by eye when 
comparing the last and penultimate AORs (again sepa- 
rated by 3.6 hr), which suggests that the charge trapping 
persistence has saturated by this time, as expected from 
the much- flattened data ramp seen in Figure 1. 

The IRAC ramp is known to exhibit a behavior which 
depends on the level of illumination, with more intensely 
illuminated pixels exhibiting a steeper initial ramp and 
saturating more quickly (these pixels' charge traps are 
filled more quickly because more free photoelectrons are 
available). We see a hint of this behavior in our data. 
Though pointing variations prevent us from tracking 
the response of individual pixels, we extract photometry 
(again via PSF fitting) using both 3- and 5-pixel-wide 
square apertures. The 3-pixel photometry - which is 
weighted somewhat more heavily by the most intensely 
illuminated pixels than is the 5-pixel photometry - shows 
a hint of a steeper ramp. We take this as further tentative 
support for our hypothesis that our ramp has a common 
origin with the IRAC ramp. The ramp behavior remains 
unchanged when we use a wider aperture, but this may 
not be diagnostic since the gradient in illumination level 
quickly flattens out beyond a few pixels. 

We would like to know why we see this ramp, espe- 
cially considering that no previous MIPS observations 
detected this effect. However, we can find no consistent 
discriminant between the presence or absence of a ramp 
in MIPS data and the state of either instrument or obser- 
vatory. The first set of AORs in ClO's observations (the 
first ~ 10 hr) were anomalously low (^ 0.3%) compared 
to subsequent observations, which they attributed to a 
thermal anneal of the 24 /xm detector conducted < 1 h 
before these observations^'^. No ramp was observed in 
the continuous, long-duration MIPS observations of ei- 
ther Knutson ct al. (2009b) or CIO, which were taken 
> 1 day after the last 24 /im anneal. The photometry 
shown in Figure 1 also occurred > 1 day after the last 
24 /im anneal, so annealing seems unlikely to explain the 
presence of the ramp in our data. 

We investigated whether preflashing could explain the 
absence of any ramp in other MIPS phase curve obser- 
vations. To preflash is to conduct a set of brief (< 1 hr) 
observations of a bright target before observing a fainter 
exoplanet system (Seager & Deming 2009; Knutson et al. 
2011); experience shows that this tends to reduce the am- 
plitude of the ramp, presumably by partially saturating 
the detector's charge traps. HD 209458 is the faintest 
of the three exoplanet systems with long-duration MIPS 
24 /im observations, but the flux difference (~20 mJy for 
HD 209458 vs. ~ 60 mJy for HD 189733) does not seem 
sufficiently large for only one of our five observations of 
HD 209458 to fail to pre-flash the detector. If the dif- 
ference were due to the increased flux from HD 189733, 
we should still see a shorter, steeper ramp at the start 

As recorded in the Spitzer observing logs, available 
at http : / /ssc . spitzer . caltech. edu/wariimiiss ion/ scheduling/ 
observinglogs/ 



of these observations. That no ramp has been reported 
previously, and that we see a ramp in the HD 209458 
data only intermittently, suggests that some other phe- 
nomenon may be at work here. 

The phase curve observations of both HD 189733 and 
HD 209458 began immediately after a data downlink to 
Earth, so this factor also does not distinguish between 
the cases. Prior to the data downlinks, our 2008 obser- 
vations of HD 209458 were preceded by 24 /im observa- 
tions of the faint RXCJ0145. 2-6033 (~ 4 mJy), but no 
24 /im observations whatsoever were made in the ~day 
leading up to Knutson ct al. (2009b) 's observations of 
HD 189733. While MIPS was operational all its ar- 
rays were continuously exposed to the sky: although 
the Spitzer operations staff planned observations so as 
to avoid placing bright sources on the 24 /im array (us- 
ing IRAS 25 /im images as a guide; A. Noriega-Crespo, 
private communication) we cannot dismiss the possibil- 
ity that occasionally some bright sources may have been 
missed. 

Thus we cannot conclusively determine why the MIPS 
observations of HD 209458 we present here show the 
detector ramp while previous, comparable observations 
have not shown such an effect. Nonetheless, the sim- 
ilarity between our photometry in Figure 1 and raw 
IRAC 8/im photometry (e.g., Agol et al. 2010) strongly 
suggests that the most likely explanation involves detec- 
tor response variations due to charge-trapping. 

3.2. The Fallback 

After the detector ramp, the photometry in Figure 1 
decreases over the rest of the observations by ~ 0.2%; we 
term this flux diminution the "fallback." The amplitude 
of this effect is of the approximate amplitude expected for 
a 24 /im planetary thermal phase curve (Showman ct al. 
2009; Burrows ct al. 2010), so our first inclination was 
to ascribe a planetary origin to this flux decrease. How- 
ever, there is a distinct qualitative similarity between 
the phase curve photometry and pre-launch calibration 
data taken with the MIPS 24 /im detector under bright 
(170 MJy sr^^) illumination, shown in Figure 2 (repro- 
duced from Young ct al. 2003). A comparison of this 
figure and Figure 1 reveals that both display the same 
qualitative signature of an early, steep ramp followed by 
a slow, gradual fallback in measured flux. The only dif- 
ferences are (1) an initial steep decrease in flux in the 
calibration data not seen in our stellar photometry (at- 
tributed by Young ct al. 2003, to the response of the 
detector to a thermal anneal immediately preceding the 
data shown) , and (2) longer ramp and fallback time con- 
stants in our data set. 

The brightest pixels in the HD 209458 MIPS obser- 
vations reach a flux of 45 MJy sr^^ (corresponding to 
1000 DN s"^). Perhaps, hke in some preflashed IRAC 
observations (cf. Knutson et al. 2011), the lower illumi- 
nation level in the HD 209458 photometry (relative to 
the stimulation response curve from Young et al. 2003) 
is responsible for the different timescales evident in the 
two 24 /im time series. However, the brightest pixels in 
the observations of CIO reached a flux of 9000 DN 
and no fallback is apparent in the continuous portion of 
those observations (though ClO's continuous photome- 
try did decrease monotonically by ~ 0.1 %, they demon- 
strated a coherent planetary phase curve in two data sets 
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Fig. 2. — Lab calibration data for the MIPS 24 fim array, taken 
from Young ct al. (2003; their Figure 7). The relevant data for 
comparison with Spitzer/MIPS observations are the gray diamonds 
labeled SUR (Sample Up the Ramp, the algorithm used to compute 
MIPS data numbers from pixel slopes). Young ct al. (2003) suggest 
that the initial sensitivity decrease (0—50 minutes) is related to de- 
tector response variations related to a thermal anneal immediately 
preceding the data; as we describe in Sec. 3, our data should not 
be affected by any anneal operations. The rest of the observations 
appear strikingly similar to our photometry of HD 209458, shown 
in Figure 1. For comparison with Figure 1 the peak pixel fluxes in 
the HD 209458 data frames are roughly 1000 DN/s. 

spanning several years: thus planetary emission, rather 
than an instrumental sensitivity variation, seems a more 
likely interpretation of their results). Similarly, no fall- 
back is seen in MIPS observations of HD 189733b (peak 
pixel flux - 1200 DN s"^ Knutson et al. 2009b) or of 
the fainter eclipsing M binary GU Boo (< 500 DN s"^; 
von Braun et al. 2008). Thus is seems possible that the 
fallback is linked to the presence of the detector ramp, 
which also appears only in our MIPS data set. 

We try a number of different functional forms to fit 
to the ~ 0.2% post-ramp fallback, which we fit simulta- 
neously with the ramp. These include a flat model (i.e., 
no decrease) , sinusoidal and Lambertian profiles with ar- 
bitrary amplitude and phase (representative of a plane- 
tary phase curve) , and a double-exponential of the form 
(1 — ae"*/'^!) X e~*/'^2, with r2 » n, motivated by the 
detector response variations seen in Figure 2. We de- 
cide which of these models is the most appropriate on 
the basis of the Bayesian Information Criterion (BIC^**). 
The model consisting of a ramp plus a decaying exponen- 
tial gives the lowest BIG: ^ 15 units lower than obtained 
with the sinusoidal or Lambertian models. Thus the data 
prefer an instrumental explanation for the low-level flux 
variations that we see. 

When using a sinusoidal or Lambertian model, the 
best-fit phase curve parameters describe a thermal phase 
variation which peaks well before secondary eclipse, sug- 
gesting a planetary hot spot eastward of the substellar 
point. Qualitatively, such a shift is consistent with obser- 
vations of both HD 189733b (Knutson et al. 2009b) and 
V And b (GIO). However, the phase offset determined by 
this fitting process is surprisingly large: 136 deg ±18 deg, 
a result which would seem to imply that the planet's 

BIC = kin N , where k is the number of parameters to 

be fit and A'^ the number of data points. A fit that gives a lower 
BIC is preferred over a fit with a higher BIC, and thus the BIC 
penalizes more complicated models. 



night side is hotter than its day side. Such a scenario 
has been predicted by some models (cf. Gho et al. 2003), 
but such a large phase offset is bigger than observed for 
either v And b or HD 189733b, and larger still when 
compared to expectations for this planet from more re- 
cent simulations (e.g., Rauschcr et al. 2008; Showman 
et al. 2009). We thus deem the phase curve fit with large 
offset to be an unlikely result, providing one more reason 
to doubt that the flux variation we see is of planetary 
origin. 

We also inject into the data a sinusoidal phase curve 
with zero phase offset and a peak-to-valley amplitude 
equal to our best-fit secondary eclipse depth results and 
repeat our analysis: in this case the best-fit sinusoidal 
and Lambertian models have a lower BIG value (by 12 
units) than the instrumental model, though the recovered 
amplitude and phase offset are still somewhat biased by 
the flux fallback. Although these results suggest that we 
are close to achieving the sensitivity required to constrain 
HD 209458b's thermal phase variations, our ignorance of 
the detailed morphology of the flux fallback prevents us 
from reaching a more quantitative conclusion. Thus, we 
can only conclude that the striking qualitative similar- 
ity between Figures 1 and 2 precludes us from making 
any deflnite claims as to the detection of planetary phase 
curve effects in our data. 

3.3. Instrument Stability 

As observed previously by GIO, the background flux of 
continuous MIPS photometry exhibits a roughly linear 
trend with time, with smaller, abrupt changes from one 
AOR to the next. The linear trend can be explained by 
a variation in the thermal zodiacal light as Spitzer's per- 
spective of HD 209458 changes with respect to the solar 
system, and GIO attribute the discontinuous, AOR-by- 
AOR background fluctuations to scattered light. What- 
ever the cause, these discontinuities are removed by the 
sky background subtraction, and do not appear to affect 
the final stellar photometry. 

During our 2008 observations we see a 0.5 /lA increase 
in the 24 ^m detector anneal current (MIPS data file key- 
word AD24ANLI), a decrease of 6 mK in the scan mirror 
temperature (keyword AGSMMTMP), and swings in the 
electronics box temperature (keyword AGEBOXTM) of 
up to 0.3 K. During sustained observations the electron- 
ics box appears to experience heating with some time 
lag, but with a much shorter cooling lag during obser- 
vational breaks to transmit data to Earth. Upon reex- 
amination of past observations, we find that these three 
parameters exhibit similar behavior during observations 
of HD 189733b (Kmitson ct al. 2009a) and of upsilon An- 
dromeda b (GIO). The MIPS optical train is cryogeni- 
cally cooled and separated from the non-cryogenic in- 
strument electronics (Heim et al. 1998), so it does not 
seem likely that the observed swings in the electronics 
box temperature should influence the photometry. Sim- 
ilarly, the anneal current and scan mirror temperature 
do not seem to correlate with either the ramp or the 
post-ramp flux decrease, so we conclude that these in- 
strumental variations do not affect our final photometry. 

4. TRANSITS 
4.1. Fitting Approach 
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We fit transits using uniform-disk and linear limb- 
darkened transit models (Mandel & Agol 2002), but 
(consistent with the results of Richardson et al. 2006) 
we find the limb-darkened model offers no improvement 
over the uniform-disk model (as determined by the BIC). 
We fit the transit data for: the time of center transit 
Tc^t, the impact parameter 6, the scaled stellar radius 
Ri,/a, the planet/star radius ratio Rp/R^, and the out- 
of-transit system flux F^,. We hold the period fixed at 
3.52474550 ± 0.00000018 d (Torres et al. 2008), which is 
a more precise determination than our observations can 
provide. To extract useful information from our half- 
transit event we always require that b and i?* /a have the 
same value, determined jointly from all our transits. We 
therefore perform one fit in which these two parameters 
are jointly fit, and a second fit in which we additionally 
fit jointly to Td and Rp/R^ across all transit events. 

We fit to the detector ramp in the 2008 transit by in- 
cluding a multiplicative factor of the form 1 — ae"*/"^, 
where t is measured from the start of the observations. 
This formulation of the ramp model is motivated by 
a physical model of the charge-trapping phenomenon 
thought to cause the IRAC 8 fim ramp (Agol et al. 2010). 
Agol et al. (2010) find a ramp based on two exponentials 
to be preferred for their high S/N observations, but we 
find that our data are not precise enough to constrain 
this more complicated model: when fitting a double- 
exponential ramp of the form 1 — aie"*/"^^ — 026"*^"^^ 
(Agol et al. 2010) the parameters for the two exponen- 
tial trends become degenerate, and the resulting fits are 
not preferred to the single ramp fit on the basis of the 
BIC. Finally, we include in all our fits the fourteen sensi- 
tivity correction terms (q) corresponding to the fourteen 
MIPS dither positions. 

4.2. Results 

Table 2 lists the results of the fit in which we assume 
a constant orbit and transit - holding b, R^,/a, Rp/R^,, 
and Tc_t constant across all transits - while Table 3 lists 
the results of the fit in which Rp/R^, and T^t (but not 
b or Ri,/a) are allowed to vary between events. We plot 
the results of fits to each individual transit, and to the 
combined data set, in Figure 3. We show how the resid- 
uals to the combined fit bin down with increasing sample 
size in Figure 4: the curve shown tracks closely with the 
N~^/'^ expectation from uncorrelated noise on short time 
scales (< 20 min), but on longer time scales the residu- 
als bin down more slowly than this. This indicates the 
presence of correlated (red) noise (cf. Pont et al. 2006) in 
these data, which is not surprising considering the ramp 
residuals apparent in Figure 3. 

We examine the residuals to the fourteen individual 
channels and see some evidence for qualitatively differ- 
ent correlated noise at different dither positions. We do 
not think it likely that this behavior is related to an 
intrapixel effect (as observed in IRAC; cf. Charbonneau 
et al. 2005), because the residual behavior we see does not 
correlate with mean PSF position relative to the bound- 
aries of individual pixels. Instead, it seems more likely 
to be a manifestation of the known position-dependent 
sensitivity effect previously attributed to residual flat- 
fielding errors (Crossfield ct al. 2010). 

The resulting posterior distributions are all unimodal 



(except for the impact parameter 6), and the usual 
correlations are apparent between h and R^,/a and be- 
tween and Rp/R^ (cf. Burke et al. 2007). As noted 
above, the 2008 transit data are strongly affected by 
the detector ramp, and we see correlations between the 
ramp parameters and the transit depth. We compute 
the two-dimensional posterior distributions of Rp/R^, 
T, and a (marginalized over all other parameters) from 
the MCMC chains using the kernel density estimate ap- 
proach described in CIO; we show these distributions in 
Figure 5 and list the elements of these parameters' co- 
variance matrix in Table 4. 

4.3. Discussion 

The three independently-fit transit depths listed in Ta- 
ble 3 have a fractional dispersion of 3%, consistent with 
our individual uncertainty estimates of 3-10%. We thus 
find no evidence for variations in transit depth, and our 
transit depths are consistent with the depth measured 
from the combination of our first two transit data sets 
(Richardson et al. 2006). 

We plot the ensemble of HD 209458b's transit depth 
measurements in Figure 6 along with a model of tran- 
sit depth vs. wavelength from Fortney et al. (2010). 
The model is consistent with the 24 /xm measurement 
we present here and agrees fairly well with the optical 
measurements of Sing ct al. (2008) and the IRAC 3.6 
and 4.5 fj,m measurements of Beaulieu ct al. (2010). How- 
ever, our model strongly disagrees with the IRAC 5.8 
and 8.0 yitm, which was also shown for the same HD 
209458b model in Fortney et al. (2010). The large dis- 
crepancy remains unclear. Given the known wavelength- 
dependent water vapor opacity, Shabram et al. (2011) 
showed that reaching all four 4 IRAC data points may 
be impossible within the framework of a simple transmis- 
sion spectrum model. Our transmission spectrum meth- 
ods are described in these papers, and the atmospheric 
pressure-temperature profile is from a planet-wide aver- 
age no-inversion model shown in Figure 7. 

We resample the posterior distributions of the indepen- 
dent transit ephemerides shown in Table 3 to determine 
our own, independent constraint on the planet's orbital 
period (assuming it is constant) using a linear relation. 
We compute the center-of-transit time and period to be 
2453549.2075 ± 0.0013 d and 3.5247537 ± 0.0000049 d, 
respectively; the covariance between these two parame- 
ters is -5.652 X 10-9 d^ The period we obtain differs 
from the established period (Torres et al. 2008) by only 
8.2 X 10"^ d (0.71 s), well within the uncertainties. 

5. SECONDARY ECLIPSES 

5.1. Fitting Approach 

We fit secondary eclipses using the uniform-disk occul- 
tation formulae of Mandel & Agol (2002), fitting each 
event for three astrophysical parameters: time of center 
of eclipse Tc^e, stellar flux F*, and eclipse depth Fp/F^ 
- as well as the fourteen sensitivity correction terms (q) 
discussed previously. We hold all other other orbital pa- 
rameters flxed at the values listed in Torres et al. (2008), 
which are more precise than our constraints based on the 
24 /zm transit photometry. We perform four different fits: 
an independent fit of each eclipse taken in isolation, and 
a fit to the combined data set in which we fit for a single 
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Fig. 3. — MIPS 24 fim transits of HD 209458b. The top panels show photometry and the best-fit model, and the lower panels show the 
residuals to the fits. For plotting purposes the data have been corrected for the MIPS 14-position sensitivity variations, normalized by the 
stellar flux, and binned by 70 points (for the individual transits) and by 210 points (for the combined data set). We also corrected for 
the ramp in the 2008 data set (corrected, binned data shown as small points) before combining the data to plot the data in the rightmost 
panel. 




1 10 100 

Bin Size [minutes] 
Fig. 4. — Dispersion of the binned residuals (solid lines) to the 
combined transit and eclipse light curve fits shown in Figure 3 
and 8. On longer timescales both fits exhibit a binned disper- 
sion 10-30% higher than expected from uncorrelated noise (dashed 
line). The dashed lines show the expectation for uncorrelated er- 
rors, which scale as N~^^'^ . The vertical dotted line indicates the 
transit duration. 



eclipse depth, but still allow Tc^e and to vary for each 
event. We use only a subset of the long-duration phase 
curve observations to fit the 2008 eclipse, as indicated in 
Table 1. We tried including a linear slope in the com- 
bined eclipse fit, but this extra parameter is not justified 
because it gives a higher BIG than fits without such a 
slope. 

5.2. Results 

The parameters for the fit in which Tc^e and Fp/F^ 
are fit jointly across all eclipses (but F^ remains inde- 
pendent) are shown in Table 5, and parameters for the 
three wholly independent eclipse fits are shown in Ta- 
ble 6. The data, best fit models, and residuals for all 
three eclipses and the combined data set are plotted in 
Figure 8. The only strong correlations apparent in the re- 
sulting posterior distributions are between F# and Fp/F^ 




Fig. 5. — Posterior distributions of the ramp parameters {a, t) 
and Rp/Rt, estimated from the MCMC analysis of the 2008 transit 
data. The ' X ' symbols indicate the best-fit parameters listed in 
Table 3, and the lines indicate the 68.27%, 95.45%, and 99.73% 
confidence intervals. The elements of these parameters' covariance 
matrix are listed in Table 4. 



Solar-abundance model 
Sing et al, 2003 
Beaulieu et al. 2010 
This work 




Fig. 6.— Measurements of 1;1icf');?an^iT'depth of HD 209458b: 
binned optical spectroscopy (Sing et al. 2008), previous mid- 
infrared photometry (Beaulieu et al. 2010), and our 24 ^m measure- 
ment. The solid line is a model generated using the (dot-dashed) 
temperature-pressure profile shown in Figure 7. The solid black 
points without errorbars represent the weighted averages of the 
model over the corresponding bandpasses (indicated at bottom). 



- expected since we are making a relative measurement. 
We show how the residuals to the combined fit bin down 
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Temperature [K] 

Fig. 7. — Temperature-pressure (T-P) profiles used to generate 
our model spectra. The dot-dashed curve is a planet-wide average 
T-P profile taken from a full (47r) redistribution model, and is used 
to model the tranmission spectrum shovirn in Figure 6. It includes 
TiO/VO opacity, but these species have only a minor effect since 
nearly all of the Ti/V has condensed out of the gas phase at these 
cooler temperatures. The solid curve is from a model assuming no 
redistribution of absorbed energy (making it hotter), and includes 
TiO/VO to drive the temperature inversion seen in Figure 9. 

with increasing sample size in Figure 4: the residuals 
average down more slowly than the N~^/'^ expectation 
from uncorrelated errors. This indicates the presence of 
correlated (red) noise (cf. Pont et al. 2006) in these data, 
which is expected given the behavior of the eclipse resid- 
uals shown in Figure 8. 

5.3. Discussion 

The three eclipse depths have a dispersion of 13%, con- 
sistent with our estimated measurement errors (12-18%). 
We thus find no evidence for variability of planetary emis- 
sion, in good agreement with general circulation models 
which predict HD 209458b's MIR dayside emission will 
vary by < 5% (e.g., Rauscher ct al. 2008; Showman ct al. 
2009; Dobbs-Dixon et al. 2010) and consistent with the 
measurement that HD 189733b's 8 /im dayside emission 
varies by < 2.7%(Agol et al. 2010). Our mean eclipse 
depth over all three epochs - 0.338% ±0.026% - is - 1.3cr 
deeper than the initial measurement by Deming et al. 
(2005) of 0.26% ± 0.046%. We convert this eclipse depth 
to a brightness temperature of 1320 K ± 80 K using the 
method outlined by CIO. 

We plot the ensemble of HD 209458b's secondary 
eclipse measurements in Figure 9 along with a model of 
planet/star contrast ratio vs. wavelength. The modeling 
procedure is described in detail in Fortney ct al. (2006) 
and Fortney et al. (2008). Using a stellar model for the 
incident flux and a solar metallicity atmosphere, we de- 
rive a radiative-convective pressure-temperature profile 
assuming chemical equilibrium mixing ratios. The model 
assumes no loss of absorbed energy to the night side, and 
redistribution of energy over the day side only (see Fort- 
ney et al. 2008). We show the pressure-temperature pro- 
file, which feature a temperature inversion due to the ab- 
sorption of stellar fiux by TiO and VO gasses, in Figure 7. 
Clearly a stronger temperature inversion is needed, as 
the contrast between the IRAC 3.6 and 4.5 /zm bands is 
not large enough. Since the 24 /im photosphere is pre- 
dicted to lie at 1 - 10 mbar on HD 209458b (Showman 
ct al. 2009) our measurement indicates a somewhat cooler 
temperature than is expected for this planet given its at- 



mospheric temperature inversion. The anomalously low 
24 /xm fiux has been noted previously (e.g., Madhusudhan 
& Seager 2010); taken in concert with v And b's large and 
still-unexplained 24 /im phase offset (CIO) these results 
suggest that our current understanding of atmospheric 
opacity sources in this wavelength range may be incom- 
plete. Alternatively, we can fit reasonably fit the 3.6, 
8.0, and 24 /im points with the dayside emission of the 
3D general circulation model of Showman et al. (2009), 
which is cooler than the corresponding ID model from 
Fortney ct al. (2008). Clearly more work is needed to 
robustly fit the dayside photometry of the planet within 
the framework of a ID or 3D self-consistent model. 

We also fit a linear relation to the three eclipse times 
in the same manner as in Sec. 4. We compute a pe- 
riod of 3.5247445 ± 0.0000097 d, which differs from the 
established period (Torres et al. 2008) by 0.99 x 10"^ d 
(0.086 s), well within the uncertainties. This value also 
agrees with our measurement of the period from the tran- 
sit fits; the two periods differ by only 0.79 ±0.94 s, which 
is (as expected) consistent with zero. 

6. JOINT ORBITAL CONSTRAINTS AND SYSTEM 

FLUX 

6.1. Timing and Eccentricity: Still a Chance for 
Winds 

Measuring the times of transit and secondary eclipse 
constrains the quantity e cos w, where e is the planet's or- 
bital eccentricity and w its longitude of periastron (Sea- 
ger 2011, chapter by J. Winn). We resample the posterior 
distributions of Tc^t and T^^e from the fits shown in Ta- 
bles 2 and 5 and compute the difference between our tran- 
sit and eclipse ephemerides (i.e., [T^^e — ?c,t] mod P—^) 
to be 32 ± 129 s after also accounting for the 47 s light 
travel time from the planet's location during eclipse to 
its location during transit (Torres ct al. 2008). This re- 
sults constrains ecosw to be 0.00004±0. 00033, consistent 
with zero and with previous constraints from radial ve- 
locity (Torres et al. 2008). We do not see the marginal 
timing offset previously reported (Knutson et al. 2008), 
which may have been biased by the higher level of corre- 
lated noise (due to the IRAC intrapixel effect) in the 3.6 
and 4.5 /im IRAC data. 

A measurement of ecosw directly constrains the ap- 
parent velocity offset that can be induced in planetary 
absorption lines during transit (cf. Montalto et al. 2011); 
this provides an independent check as to whether the re- 
cent measurement of a velocity offset of 2 ± 1 km s~^ in 
HD 209458b (Snellen et al. 2010) can be attributed to a 
low, but nonzero, orbital eccentricity. Our timing mea- 
surements of HD 209458b set a 3cr upper limit on any 
velocity offset due to the planet's orbital eccentricity of 
only 140 m / ^/Y^^~f? . Thus the claimed velocity off- 
set, though still of low significance, cannot be dismissed 
as resulting from the HD 209458b's orbital eccentricity. 

6.2. System Flux: No Excess Detected 

Although our primary science results - the transit and 
eclipse depths - rely on relative fiux measurements, our 
observations also allow us to measure absolute 24 //m 
photometry for the HD 209458 system. Our fiux mea- 
surements for this system vary from epoch to epoch 
by much more than our quoted statistical uncertainties. 
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Fig. 8. — MIPS 24 /im eclipses of HD 209458b. The top panels show photometry and the best-fit models, and the lower panels show the 
residuals to the fits. For plotting purposes the data have been corrected for the MIPS 14-position sensitivity variations, normalized by the 
stellar fiux, and binned by 70 points (for the individual eclipses) and by 210 points (for the combined data set). 

We have described a homogeneous analysis of 
all Spitzer MIPS observations of the hot Jupiter 
HD 209458b. The data comprise three eclipses, two and 
a half transits, and a long, continuous observation de- 
signed to observe the planet's thermal phase curve; of 
these, analysis of two of the eclipses, one transit, and 
the phase curve observations have remained unpublished 
until now. The long-duration phase curve observation ex- 
hibits a detector ramp that appears similar to the ramp 
seen in Spitzer/IRAC 8 /im photometry, and we model 
this effect using the exponential function proposed by 
Agol et al. (2010). We also see a ~ —0.2% flux decrease in 
the latter portion of the phase curve observations. This 
fallback is similar to a known (but poorly characterized) 
variation in the response of the MIPS detector when sub- 
jected to bright illumination (cf. Figure 2 and Young 
et al. 2003). 

We are unable to determine why either the fallback or 
the ramp have not been seen in any prior MIPS observa- 
tions. Despite this failure the correspondence between 
our photometry and the pre-launch array calibration 
data leads us to conclude that ramp and fallback are cor- 
related and both are most likely of instrumental, rather 
than astrophysical, origin. This conclusion is strength- 
ened by the result that fitting periodic phase functions 
to the data yields a planetary model hotter on its night 
side than its day side, strikingly at odds with theory 
(Rauscher et al. 2008; Showman et al. 2009; Dobbs-Dixon 
ct al. 2010; Cowan & Agol 2011a) and inconsistent with 
other published observations of hot Jupiters (Harrington 
ct al. 2006; Knutson et al. 2007, 2009b; Crossfield ct al. 
2010). 

We see no evidence for variation in the three eclipse 
depths, and a joint fit of all three eclipses gives our best 
estimate of the 24 /im planet /star contrast: 0.338% ± 
0.026%. This value is more precise and higher than the 
previously published measurement (Doming et al. 2005), 
and corresponds to an average dayside brightness tem- 
perature of 1320 K±80 K, consistent with models of this 
planet's thermal atmospheric structure (Showman et al. 
2009; Madhusudhan & Seager 2010; Moses et al. 2011). 
We note parenthetically that this new eclipse depth has 
already diffused into several papers (cf. Showman et al. 
2009; Burrows ct al. 2010; Madhusudhan & Seager 2009, 
2010; Fortney et al. 2010; Moses et al. 2011); the value 



Fig. 9.— Measurements W^fiPia'^^fllldary eclipse depth of 
HD 209458b: previous Spitzer/IRAC photometry (Knutson ct al. 
2008) and our 24 /^m measurement. The solid line is from a model 
assuming zero redistribution of incident flux and including gaseous 
TiO and VO to drive a temperature inversion; we show this model's 
temperature-pressure proflle in Figure 7. The dashed line is the 
emission spectrum from Showman ct al. (2009). The solid black 
points without errorbars represent the weighted averages of the 
models over the corresponding bandpasses (indicated at bottom). 



but the variations are not large compared to the < 1% 
repeatability and 2% absolute calibration accuracy of 
the MIPS 24 /zm array (Engclbracht et al. 2007). Our 
21 pixel aperture encloses 99.2% of the starlight (as de- 
termined from our synthetic PSF), and we account for 
this small effect in the value quoted below. 

We therefore report the 24 /im system flux as 18.7 ± 
0.4 mjy, consistent with the flux expected from the 
HD 209458 stellar photosphere (as reported by Dom- 
ing et al. 2005). HD 209458 was not detected by IRAS 
(Beichman ot al. 1988), but is present in the Wideflold 
Infrared Survey Explorer's all-sky point source catalogue 
(Wright ot al. 2010). The WISE photometry gives a W4 
system flux of 25.74±0.12 mJy, which is higher than but 
marginally 3ct) consistent with the Spitzer-derived 
value after accounting for the different wavelengths of the 
two instruments. We therefore conclude that HD 209458 
does not have a strong 24 /im infrared excess, as is typical 
of middle-aged F dwarfs (Moor ct al. 2011). 

7. CONCLUSIONS AND FUTURE WORK 
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and uncertainty quoted in those works are close to those 
we report here, so their conclusions should be relatively 
unaffected. 

We see no evidence for variations in our transit mea- 
surements, and a joint fit of our two and a half tran- 
sits yields a 24 ^m transit depth of 1.484% ± 0.033%. 
The transit depth is less well-constrained than the eclipse 
depth because only half of the first transit was observed, 
and the last transit occurred during the detector ramp. 

The ephemerides calculated from our analyses of the 
transits and eclipses allow us to compute orbital periods 
of 3.5247537 ±0.0000049 d and 3.5247445 ±0.0000097 d, 
respectively, which are consistent with but less precise 
than the orbital period of Torres ct al. (2008). Eclipses 
occur 32±129 s earlier than would be expected from a cir- 
cular orbit, which constrains the orbital quantity ecosw 
to be 0.00004±0.00033. This suggests that HD 209458b's 
inflated radius (larger than predicted by models of plane- 
tary interiors; Fortney et al. 2007) cannot be explained by 
interior heating from ongoing tidal circularization, and 
that the possible velocity offset reported by Snellen ct al. 
(2010) cannot be explained by a nonzero orbital eccen- 
tricity. 

Although we obtain improved estimates of the 24 fim 
transit and secondary eclipse parameters, instrumental 
effects prevent a conclusive detection of the planet's ther- 
mal phase curve. The phase curve signal is inextricably 
combined with the systematic fallback effect, despite es- 
timates that the planet's day/night contrast should be as 
large as a few parts per thousand (Showman et al. 2009). 
Such a large and intermittent systematic effect has pro- 
found implications for future mid-infrared exoplanet ob- 
servations with EChO, SPICA, and JWST. Models of 
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terrestrial planet phase curves predict phase amplitudes 
of< lO""* (Selsiset al. 2011; Maurin et al. 2011); such ob- 
servations could be utterly confounded by the instrumen- 
tal systeniatics seen in our observations, and so may be 
much more challenging that has been heretofore assumed 
(Kahcnegger & Traub 2009; Seager & Deming 2009). Al- 
though it may be possible to reduce the effect of the ramp 
with a pre-flash strategy similar to that adopted for the 
8 /im IRAC array, a further defense against these chal- 
lenges would seem to be a more comprehensive campaign 
of array characterization. Specifically, a detailed char- 
acterization of the detector response to sustained levels 
of the high illumination expected from observations of 
terrestrial planets around the brightest nearby stars is 
highly desirable, and should be considered an essential 
requirement for all future infrared space telescopes. 
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TABLE 1 

Spitzer /MIPS 24 /iM Observations of HD 209458b 



UT Date 



Event Duration (hr) ti„t (s) N 



exposures 



Bkd (MJy Sr-i)'' At (s) 



UT 2004 Dec 5 
UT 2004 Dec 6 
UT 2005 Jun 27 
UT 2005 Dec 1 
UT 2008 Jul 25= 
UT 2008 Jul 27" 



Half transit 
Eclipse 
Transit 
Eclipse 
Transifl 
Eclipse 



2.8 
5.8 
5.6 
5.6 
14.2 
6.9 



8.91 
9.96 
9.96 
9.96 
9.96 
9.96 



840 
1680 
1680 
1680 
4060 
2072 



28.8 
29.2 
28.9 
26.7 
28.3 
27.9 



-544 
-531 
-612 
+183 
-649 
-663 



^ Average sky backgrounds as reported by DRIBKGND keyword. 
^ For each event, At = (HJD) - {B.JDtdb) 

" These events were observed as part of a single, continuous phase curve observation with a duration of 58 
hours spanning one transit and one secondary eclipses. 

^ This transit was corniplcd by an apparent ramp in detector sensitivity, so we used a longer section of data 



to better constrain the ramp parameters in the joint fit. 



TABLE 2 
Joint Transit Fits 



Parameter 



2004 



2005 



2008 



CO 

ci 

C2 
C3 
C4 
C5 
C6 
C7 
C8 

eg 

CIO 
Cll 
C12 
CI 3 

F* [mjy] 
a 

T[d] 



+0.0009 ± 0.0022 
+0.0097 ± 0.0027 
-0.0030 ± 0.0022 
+0.0106 ± 0.0024 
-0.0034 ± 0.0024 
-0.0021 ± 0.0032 
+0.0017 ± 0.0022 
-0.0023 ± 0.0033 
-0.0036 ± 0.0023 
+0.0098 ± 0.0023 
-0.0032 ± 0.0028 
-0.0126 ± 0.0032 
+0.0008 ± 0.0028 
-0.0028 ± 0.0024 
+18.845 ± 0.012 



+0.00036 ± 0.00086 
+0.00881 ± 0.00088 
+0.00101 ± 0.00136 
+0.00866 ± 0.00079 
+0.00148 ± 0.00115 
+0.01006 ± 0.00080 
-0.00574 ± 0.00078 
+0.00167 ± 0.00081 
-0.00298 ± 0.00078 
-0.00324 ± 0.00135 
+0.00394 ± 0.00111 
-0.00974 ± 0.00103 
-0.00039 ± 0.00085 
-0.01360 ± 0.00091 
+18.7784 ± 0.0049 



+0.00274 ± 0.00053 
+0.01387 ± 0.00055 
+0.00487 ± 0.00053 
+0.00858 ± 0.00056 
-0.00020 ± 0.00053 
+0.01310 ± 0.00052 
-0.00502 ± 0.00057 
-0.00720 ± 0.00054 
-0.00121 ± 0.00058 
-0.00919 ± 0.00060 
+0.00210 ± 0.00057 
-0.01016 ± 0.00057 
-0.00122 ± 0.00065 
-0.01062 ± 0.00052 

+18.696 ± 0.010 
+0.02437 ± 0.00068 

+0.174 ± 0.016 



[BJDtdbJ 
b 

Rt/a 
(Rp/R*)'^^ 



2453549.20852 ± 0.00049 
+0.590 ± 0.062 
+0.1205 ± 0.0066 
+0.1218 ± 0.0014 
+0.01483 ± 0.00033 



Computed from the posterior MCMC distributions of Rp /R* . 
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TABLE 3 
Semi-Joint Transit Fits 



Parameter 


2004 


2005 


2008 


Co 


-U.UUUi ± U.UUz4 


+U.UUU40 ± U.UUUo4 


+U.UUZDO ± U.UUUo4 


Cl 


+0.0091 ± 0.0030 


-1-0.00866 ± 0.00095 


+0.01368 ± 0.00053 


C2 


n r\n'DO i n r\noo 

-U.UUoi ± U.UUzJ 


H-U.UUUo7 ± U.UU149 


+(J.UU524 ± U.UUUyd 


C3 




H-U. UUo/z it U.UUUoZ 


H-U.UUoro it U.UUUuu 


C4 


-U.UUoz it U.UUzo 


4-U.UU14D it U.UUilU 


-U.UUUio it U.UUUoo 


C5 


n 0091 -1- n 0090 
-u.uuzi m u.uuzy 


J_0 0100Q -1- 0007)>i 


4-0 01 "^1 4- OOOf^*? 


C6 


4-0 001 Q -1- 0099 


OOfiO*^ A- OOOkil 


OO'^I fi 4- 000'^'^ 


c? 


n 001 Q -1- 00*^0 


4-0 001 '^y 4- 0008*^ 


OOfif^fi 4- OOO'^J^ 


Cs 


-0.0039 ± 0.0023 


-0.00301 ± 0.00078 


-0.00110 ± 0.00063 


C9 


-1-0.0095 ± 0.0022 


-0.00317 ± 0.00144 


-0.00910 ± 0.00057 


ClO 


-0.0032 ± 0.0029 


-f0.00378 ± 0.00093 


+0.00189 ± 0.00052 


Cll 


-0.0124 it 0.0028 


-0.00942 it 0.00082 


-0.00977 ± 0.00056 


C12 


-1-0.0007 ± 0.0028 


-0.00035 it 0.00086 


-0.00189 ± 0.00129 


Cl3 


-0.0026 ± 0.0022 


-0.01330 it 0.00078 


-0.01096 ± 0.00061 


Tc,t [BJDtds] 


2453344.7718 ± 0.0025 


2453549.20746 it 0.00065 


2454673.60391 ± 0.00074 




-1-0.1227 ± 0.0060 


-to. 1189 it 0.0020 


+0.1238 ± 0.0019 


F* [mjy] 


-H8.850 ± 0.016 


+18.7735 ± 0.0056 


+18.6947 ± 0.0095 




-1-0.0151 ± 0.0015 


+0.01413 ± 0.00046 


+0.01531 ± 0.00046 








+0.02458 ± 0.00071 


r [d] 






+0.168 ± 0.016 




b 

R,/a 


+0.581 ± 0.070 
+0.1197 ± 0.0069 





^ Computed from the posterior MCMC distributions of Rp/R^. 



TABLE 4 
Ramp and transit 
depth covariance 
MATRIX (CF. Figure 5). 



Element Value/10^ 



0.507 

-'[d] 266 d2 

3-49 

o"a,T[d] -5.03 d 

'^a,Rp/R, 0-379 

'^T[ti],i?p/B. -16.3 d 



TABLE 5 
Joint Eclipse Fits 



Parameter 


2004 


2005 


2008 


CO 


+0.00298 ± 0.00080 


+0.00048 ± 0.00078 


+0.00278 ± 0.00073 


Cl 


+0.00793 ± 0.00078 


+0.01103 ± 0.00080 


+0.01331 ± 0.00071 


C2 


+0.00619 ± 0.00079 


+0.00349 ± 0.00078 


+0.00397 ± 0.00074 


C3 


+0.00943 ± 0.00075 


+0.00879 ± 0.00078 


+0.00887 ± 0.00071 


C4 


+0.00048 ± 0.00108 


+0.00186 ± 0.00102 


-0.00002 ± 0.00074 


C5 


+0.01148 ± 0.00080 


+0.01151 ± 0.00079 


+0.01480 ± 0.00069 


C6 


-0.00557 ± 0.00079 


-0.00461 ± 0.00084 


-0.00432 ± 0.00072 


C7 


-0.00571 ± 0.00082 


-0.00582 ± 0.00079 


-0.00514 ± 0.00072 


C8 


-0.00171 ± 0.00078 


-0.00020 ± 0.00082 


-0.00329 ± 0.00070 


eg 


-0.00231 ± 0.00079 


-0.00361 ± 0.00079 


-0.00754 ± 0.00070 


CIO 


+0.00052 ± 0.00082 


+0.00084 ± 0.00093 


-0.00137 ± 0.00080 


Cll 


-0.00691 ± 0.00081 


-0.00973 ± 0.00076 


-0.00712 ± 0.00073 


Cl2 


-0.00507 ± 0.00082 


-0.00402 ± 0.00078 


-0.00333 ± 0.00072 


Cl3 


-0.01142 ± 0.00079 


-0.00968 ± 0.00081 


-0.01121 ± 0.00071 


F* [mJy] 


+18.78683 ± 0.00486 


+18.70529 ± 0.00485 


+18.60906 ± 0.00472 




+0.00338 ± 0.00026 






Tb [K] 


+1310 ± 80 






Fp [^Jy] 


+63.2 ± 4.9 






Te,, [BJD tdb]^ 


2453706.0595 ± 0.0014 





^ Computed from the posterior MCMC distributions of F« and Fp/F^. 

^ Jointly-fit ephemeris, assuming a period of 3.5247455 d (Torres et al. 2008). 
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TABLE 6 



Independent Eclipse Fits 


Parameter 


2004 


2005 


2008 


CO 


0.00262 ± 0.00090 


-0.00020 ± 0.00112 


0.00044 it 0.00264 


Cl 


0.00744 ± 0.00089 


0.01028 ± 0.00124 


0.01246 it 0.00114 


C2 


0.00500 ± 0.00169 


0.00396 ± 0.00084 


0.00449 it 0.00106 


C3 


0.00879 ± 0.00101 


0.00879 ± 0.00077 


0.00850 it 0.00075 


C4 


-0.00042 ± 0.00090 


0.00198 ± 0.00111 


0.00143 it 0.00160 


C5 


0.01090 it 0.00104 


0.01109 it 0.00091 


0.01392 it 0.00115 


C6 


-0.00488 ± O.OOlOo 


-0.00502 ± 0.00109 


-0.00351 ± 0.00120 


C7 


-0.00453 ± 0.00138 


-0.00331 it 0.00284 


-0.00561 it 0.00086 


eg 


-0.00248 ± 0.00109 


-0.00078 it 0.00089 


-0.00249 it 0.00105 


eg 


-0.00230 ± 0.00082 


-0.00239 it 0.00142 


-0.00744 it 0.00073 


CIO 


0.00113 ± 0.00118 


0.00086 it 0.00093 


0.00079 it 0.00225 


Cll 


-0.00588 ± 0.00118 


-0.01021 it 0.00097 


-0.01012 it 0.00328 


Cl2 


-0.00315 ± 0.00201 


-0.00385 it 0.00078 


-0.00189 it 0.00171 


C13 


-0.01197 ± 0.00104 


-0.01090 it 0.00166 


-0.01059 it 0.00098 


Tc,e [BJD tdb]^ 


2453346.5348 ± 0.0028 


2453706.0600 it 0.0029 


2454675.3639 it 0.0026 




0.00325 ± 0.00053 


0.00384 it 0.00046 


0.00281 it 0.00051 


Tb [K] 


1270 ± 190 


1450it 230 


1130 it 160 


F, [mjy] 


18.7884 ± 0.0069 


18.7023 it 0.0064 


18.6155 it 0.0068 


Fp [/^Jy] 


61.1 ± 10.0 


71.9 it 8.6 


52.3 it 9.5 



^ Jointly-fit cphcmcris, assuming a period of 3.5247455 d (Torres et al. 2008). 
^ Computed from the posterior MCMC distributions of i^* and Fp/F^. 



